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We report on a search for CP violation in r — > Kgir v T decays using a data sample of 699 fb -1 
collected by the Belle experiment at the KEKB electron-positron asymmetric-energy collider. The 
CP asymmetry is measured in four bins of the invariant mass of the system and found to 

be compatible with zero with a precision of O(10 -3 ) in each mass bin. Limits for the CP violation 
parameter Q(ns) are given at the 90% confidence level. These limits are |Sj(?7s)| < 0.026 or better, 
depending on the parameterization used to describe the hadronic form factors, and improve upon 
previous limits by one order of magnitude. 

PACS numbers: 13.35. Dx, 11.30.Er, 14.80. Fd 



To date CP violation (CPV) has been observed only 
in the K and B meson systems. In the Standard Model 
(SM), all observed CPV effects can be explained by the 
irreducible complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix To find new 
physics, it is important to look for other CP-violating ef- 
fects in as many systems as possible. One such system 
is the r lepton. In hadronic r decays, one can search for 
CPV effects of possible new physics that could originate, 
for example^ from the Minimal Supcrsymmetric Standard 
Model 0, y] or from multi-Higgs-doublet models 0, [H[ 
that play an important role in strangeness changing pro- 
cesses. 

This paper describes a search for CPV in t ± — s- 
Kgir^v-r decays. It should be noted that CPV in K° 
decays leads to a small SM CP asymmetry of O(10~ 3 ) in 
the rates of this r decay mode [|| . This asymmetry is 
just below our experimental sensitivity. Here the focus 
will be on CPV that could arise from a charged scalar bo- 
son exchange @], e.g., a charged Higgs boson. This type 
of CPV cannot be observed from measurement of ^ de- 
cay rates. However, it can be detected as a difference 
in the decay angular distributions and is accessible 
without requiring information about the t polarization 
or the determination of the r rest frame. Limits for the 
CPV parameter in this decay mode have been published 
previously by the CLEO collaboration from an analysis 
of 13.3 fb" 1 of data @. 

In the SM, the differential decay width in the hadronic 
rest frame (<fi + q? = 0) is given by (see Q for details) 
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i K - Q 2 f 



(47T) 



I9i I 



(1) 



dQ 2 dcos 9 da dcos (3 
J755 2 2^ 2 ' 



where Gf is the Fermi coupling constant, 9 C is the 
Cabibbo angle, m T is the mass of the r lepton, q\ and 
g*2 denote the three-momenta of the Kg and n~ , respec- 
tively, and Q 2 = (qi + q-i) 2 is the square of the invariant 
mass of the K^n^ system. The four hadronic functions 
W x with X e (B, SA, SF, SG) (see ^) are formed from 
the vector and scalar form factors F(Q 2 ) and Fs(Q 2 ) and 
are proportional to |F| 2 , |F S | 2 , 5R(FF S ), and3(FF s ), re- 
spectively. The Lx functions, which contain the angular 
dependence, can be calculated from electroweak theory 
(see [H). The angle /3 is defined by cos/3 = nr, ■ q\ where 
qi = qi/\q\ \ is the direction of the Kg and ??l is the di- 
rection of the e + e~~ center of mass (CM) system, both 
observed in the hadronic rest frame. The azimuthal an- 
gle a is not observable in this experiment and has to be 
integrated over. The variable 9 is the angle between the 
direction opposite to the direction of the CM system and 
the direction of the hadronic system in the r rest frame. 
In this experiment, the direction of the r is not known 
but 9 can be calculated from the hadronic energy 
measured in the CM system: 

2xm 2 T — m 2 — Q 2 E h 

cos9 = =, x — 2—=, (2) 

[m 2 - Q 2 )yJ\ - Am 2 T /s V s 

where s = 4E^ e denotes the squared CM energy. 

The effect of the exchange of a charged scalar boson 
can be introduced by replacing the scalar form factor F$ 
with 

F S (Q 2 ) Fs{Q 2 ) - F S (Q 2 ) + ^F H (Q 2 ), (3) 

where Fh denotes the form factor for the scalar bo- 
son exchange [Fh = (K° (qi)jT~ (q2)\us\0)} and 775 is the 
corresponding dimensionless complex coupling constant 
[1 EH ELI- The differential decay width for the CP conju- 
gate process, dT T +, is obtained from Eq. ([T]) and Eq. ([3]) 
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by the replacement r)s — > rjs- Using this relation the CP 
violating quantity is given by 0] 



A 



LW 



X 



(4) 



= |<fi|3(FF^)9(?7s) cos if) cos f3, 



where ip denotes the angle between the direction of the 
CM frame and the direction of the r as seen from the 
hadronic rest frame and can be calculated as 



cos ip 



x{m 2 T + Q 2 ) - 2Q 2 



(m 2 ~Q 2 )y/x 2 -4Q 2 /s 



(5) 



Since the CP violating term is proportional to cos j3 cos ip, 
it cancels out if one integrates over the angles (3 and 
ip, e.g., for branching fractions. Furthermore, the CP 
violating effect is only observable if ^s(FF^j) ^ 0. The 
form factor Fjj is related to the SM weak scalar form 
factor Fs via: 



Fh(Q 2 ) 



Q 2 



-F S (Q 2 ) 



(0) 



where m u and m s denote the up and strange quark 
masses, respectively. The derivation of Eq. ^ is dis- 
cussed in Q although Fh is not used there explicitly. 
The chosen value {m u — m s ) = —0.1 GeV/c 2 defines the 
scale of the CPV parameter 9(?7s). Because the CLEO 
collaboration used a different relation Fjj = AIFs with 
M = lGeV/c 2 as well as a different normalization of 
Fs(Q 2 ), 3(??s) is not the same as the CP parameter A 
that was used in In the following, the approximate 
relation 3(?7s) — — 1.1A is used to enable a comparison 
of the results. 

To extract the CP violating term in Eq. @, we de- 
fine an asymmetry in bin i of Q 2 using the difference 
of the differential r+ and r _ decay widths weighted by 
cos /3 cos ip: 



JJJn2 Z cos f3 cos if) 



>Ql 



duj 
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(cos f3 cos ip) l T -t 
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with duj = dQ 2 dcos ddcos f3. In other words, A cp is the 
difference between the mean values of cos f3 cos ip for t + 
and t~ events evaluated in bins of Q 2 . 

We use 699 ftT 1 of data collected at the T(35), T(45) 
and T(5S) resonances and off-resonance with the Belle 
detector [13j at the KEKB asymmetric-energy e + e~ col- 
lider jl4j|. The signal and backgrounds from t + t~ events 
are generated by KKMC/TAUOL A El . The detector 
response is simulated by a GEANT3 [l_6| based program. 

Using standard event topology requirements, a 



t + t (7) sample is selected as described in 17 



In the CM frame, the event is divided into two hemi- 
spheres using the plan e perpendicular to the direction 
of the thrust axis [18|. Events with one charged track 
from an electron, muon or pion in one hemisphere (tag 
side) and a charged pion and a Kg — > 7r + 7r~ candidate 
in the other hemisphere (signal side) are chosen. The 
Kg candidates are required to have an invariant mass in 
the range 0.485 GeV/c 2 < < 0.511 GeV/c 2 and a 

reconstructed K s decay length greater than 2 cm. The 
selection criteria for the signal side and particle iden- 
tification criteria are described in detail in [l9j]. Back- 
grounds from decays with a it are suppressed by reject- 
ing events containing photons on the signal side with en- 
ergies greater than 0.15 GeV. To further suppress back- 
ground from e + e _ — > qq (q — u, d, s, and c) processes, a 
thrust value above 0.9 is required and for events with a 
pion on the tag side, the number of tag side photons 
with energies greater than 0.1 GeV must be less than 
five. In total, (162.2 ± 0.4) x 10 3 t+ -> K s tt+D t and 
(162.0 ± 0.4) x 10 3 t~ -> K a s -K~v T candidates are se- 
lected. Background contributions from r decays with 
the exception of r — > v t h ± 'k + ii~ and contributions 
from e + e~ — > qq and two-photon processes are estimated 



from Monte Carlo (MC) simulation H^-Ell using the 
branching fractions from (23l |. Contributions from — > 
j/ r 7r ± 7r + 7r _ are estimated using the data in the two Kg 
sideband regions, 0.469 GeV/c 2 < < 0.482 GeV/c 2 

and 0.514 GeV/c 2 < M„ < 0.527 GeV/c 2 Q. 

The largest background contribution is due to other r 
decays, namely (9.5 ± 3.2)% of the events in the selected 
signal sample from r ± — > ^ 7 -A'gAT° 7r± 7 (3.7 ± 1.2)% from 
t± -> VrKg^ir , (1.7 ±0.2)% from r± -> iy r K s K ± , 
and (1.79 ± 0.03)% from r± -> i/^tt+tt" . The con- 
tribution from e + e _ — > qq is (3.4 ± 1.0)%. The back- 
grounds from bb, Bhabha and two-photon processes are 
negligible. The total contribution of background pro- 
cesses is (22.1 ± 3.6)%. The invariant mass of the K s tt 
system, W = y/Q 2 , for the selected data events is shown 
in Fig. Q] together with simulated signal events and the 
background contributions discussed above. Signal events 
were generated by a modified version of TAUOLA that 
incorporates the results of [l9j . 

To avoid possible bias, the CPV search is performed 
as a blind analysis. First, possible sources of artificial 
CPV, such as forward-backward (FB) asymmetries in the 
e + e~ — > t + t~ production (7 — Z interference effects and 
higher-order QED effects) and detector induced differ- 
ences between 7r + and 7r~ reconstruction efficiencies, are 
studied using data. Other unknown sources are inves- 
tigated in data by measuring the CP asymmetry in a 
control sample described below. 

The FB asymmetry is measured in — > Vt-tt^tt+tt^ 
events (excluding K s — > tt + it~ signal candidates by us- 
ing a mass and decay length veto) as a function of the 
momentum and polar angle of the tt^it+it - system. An 
effect of a few percent is observed, which is described well 
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data 

| x -» v t K s n 

t -> v, i K s K L nit" (n>0) 
| x -> v t K g ji rm° (n>0) 
I t -> v, K s K rm° 
I t -> v, it K s K s nn° (n>0) 
| other x decays 
| K s sideband data 
I no. dd, ss 



1.4 1.6 

W (GeV/c 2 ) 

FIG. 1. Mass spectrum of the Kgir^ system. Data are indi- 
cated by the squares, simulated signal and the estimated back- 
ground contributions are shown by the colored histograms. 
All background modes have been determined from Monte 
Carlo with exception of r — > u t tt ± ty + tt~ which has been 
estimated from K% sideband data. 



by the MC simulation. The asymmetry for detection, 
which can arise because of the different nuclear interac- 
tion cross sections for positively and negatively charged 
hadrons, is studied in the laboratory system as a func- 
tion of momentum and polar angle of the charged pions in 
t~^~ — > ^ T 7r ± 7r + 7r _ (excluding 7r + 7r~ combinations consis- 
tent with Kg decays) events and found to be of O(10~ 2 ) 
(see [22| for details). Using these measurements, correc- 
tion tables are obtained that are then applied as weights 
for each event. Since the CP asymmetry is measured as a 
function of angles relative to the r direction rather than 
polar angles in the laboratory, the net effect of these cor- 
rections on the CP asymmetry is very small [O(10 -4 ) for 
FB asymmetry effects and O(10~ 3 ) for the 7r ± detection 
asymmetry] . 

A control sample is selected from t ± — > i/ T 7r ± 7r + 7r _ 
events (25j by requiring that the invariant mass of both 
7r + 7r~ combinations lie outside of the Kg mass window 
but the mass of one of the combinations lie in the side- 
band of this window. The resulting sample consists of 
about 10 6 events, i.e. about three times more than the 
signal sample. The CP asymmetry measured in this con- 
trol sample is very small [O(10 -3 )] (see 27| for details) 
and serves as an estimate of the remaining unknown sys- 
tematic effects. 



The observed CP asymmetry in the selected t^- —> 
candidate sample is shown in Table U for four 
bins of the hadronic mass W = yQ^ before and after 
applying the corrections for higher-order QED and tt^ 
detection asymmetry effects. The 4th column shows the 
final values of the CP asymmetry after subtraction of the 
background contributions. Here, we assume that there 
is no CP asymmetry in the background and correct the 



TABLE I. CP asymmetry A cp measured in bins of the 
hadronic mass W. The 2nd and 3rd column show the ob- 
served asymmetry with statistical errors only, before and after 
correcting for higher-order QED and 7T detection asymme- 
try effects. The final CP asymmetry after background sub- 
traction is shown in the 4th column where first and second 
errors correspond to statistical and systematic errors, respec- 
tively. The 5th column shows the observed number of signal 
events m per W bin (after background subtraction) divided 
by N s = J2i n i- 



W 

(GeV/c 2 ) 



Observed 



A cp (10" 
Corrected 



Backgr. subtr. «i/iV s (%) 



0.625-0.890 -0.1 ±2.1 5.2 ± 2.1 

0.890-1.110 -2.7 ±1.7 1.6 ±1.7 

1.110-1.420 -5.1 ±4.7 -3.5 ±4.7 

1.420-1.775 9.3 ±12.1 9.6 ± 12.1 



7.9 ±3.0 ± 2.8 36.53 ±0.14 

1.8 ±2.1 ± 1.4 57.85 ±0.15 

-4.6 ±7.2 ± 1.7 4.87 ±0.04 

-2.3 ± 19.1 ± 5.5 0.75 ±0.02 



TABLE II. Systematic uncertainties in the CP asymmetry 
A cp . The second column shows the uncertainties due to ef- 
fects introduced by the detector, which are estimated from 
the A cp measurement in the control sample. Contributions 
from uncertainties in the background estimates and limited 
MC statistics are small in comparison. 



ir 



Systematic uncertainties (10 



(GeV/c 2 ) 


Detector 


Backgr. 


MC stat. 


Total 


0.625-0.890 


2.76 


0.59 


0.15 


2.83 


0.890-1.110 


1.40 


0.04 


0.10 


1.40 


1.110-1.420 


1.50 


0.25 


0.79 


1.71 


1.420-1.775 


5.18 


0.96 


1.38 


5.45 



background effects as 



_ (cos /3 cos tp)]^ 



A7 = 



W 6 



(cos /? cos ipy _, 



(8) 



where are the fractions of background in the selected 
r samples in W bin i. 

In order to account for possible systematic uncertain- 
ties due to detector effects, the quadratic sum of the val- 
ues of A cp measured in the control sample and their sta- 
tistical errors are used as an estimate of the systematic 
error. Other contributions to the systematic error arise 
in the background subtraction because of uncertainties 
in the estimated number of background candidates and 
limited MC statistics. These contributions are however 
small in comparison. A summary of the systematic un- 
certainties is given in Table [TXJ 

The background subtracted asymmetry is shown in 
Fig. [2] (a) and (b) with statistical and systematic errors 
added in quadrature. The asymmetry is small and except 
for the lowest mass bin within one standard deviation (cr) 
of zero. For comparison the predicted CP asymmetry is 
shown in Fig. [2] (a) for 3(?7 S ) = 0.1 and 3£(?7s) = [Ill- 
Note that the current best limit by the CLEO experiment 
Q corresponds to |5(?7s)| < 0.19. 

From the measured values of A cp the CPV parameter 
3(775) can be extracted, which allows an interpretation 
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FIG. 2. (a) Measured CP violation asymmetry after back- 
ground subtraction (squares). The vertical error bars are 
the statistical error and systematic errors added in quadra- 
ture. The CP asymmetry measured in the control sample is 
indicated by the blue triangles (statistical errors only) and 
the inverted red triangles show the expected asymmetry for 
Sr(?7s) =0.1 [3i(r?s) = 0]. (b) Expanded view (the vertical 
scale is reduced by a factor of five) . 



in the context of NP models. Taking into account the 
detector efficiencies, the relation between A cp and $s(tjs) 
is given as 



ill 

rii 



C(Q 2 



MFF* H ) 



dQ 2 = cMvs) 



(9) 



v r events 
Y]. m is the 
events. The 



where rii is the observed number of r 
in Q 2 bin i (Q 2 e [Q\^Q%^) and N s 
total number of observed r —> Ka 
function C(Q 2 ) includes the detector efficiency as well 
as all mo del- independent terms. First, the efficiency is 
determined as a function of Q 2 , f3 and 8, then C(Q 2 ) 
is obtained after numerical integration over the decay 
an gles (3 and 8. The parameterization of C(Q 2 ) is given 
in |22| . 

Using the function C(Q 2 ) and the fractions N s /rii 
which are given in Table HI the linearity constants Cj, 
which relate ^4 cp and 3(775), can be determined for any 
parameterization of the form factors F and Fg simply 
by calculating the integral in Eq. © [28j . 

To determine limits for |$J(?7s)|, three paramcteriza- 
tions of F and Fs [exploiting Eq. (j6])] as linear com- 
binations of Brcit-Wigncr shapes of the vector reso- 
nances if* (892) and i^*(1410) and the scalar resonances 
Kq(800) and 1^(1430) are used. These parameteriza- 
tions were determined in an earlier Belle measurement 
of the KgTr ± mass spectrum In addition, a con- 

stant strong interaction phase difference between F and 
F S , cj> 3 = arg[F s (Q^ in )] - arg[-F(Q^)] with Q 2 min = 
(inn + m K o ) 2 , is introduced for generality because such 
a relative phase cannot be determined from the 
mass spectrum. 

Using Eq. ([5]), the linearity constants Cj are calcu- 
lated in each mass bin for 0s = 0°, 5°, . . . , 360° and 



the obtained values of 3(7/s) with associated uncertain- 
ties are combined to determine upper limits for 
For each parameterization, the value 0s giving the most 
conservative limit is chosen. For the three parameter- 
izations of F and Fs, this results in the range of lim- 
its \3(ris)\ < (0.012 - 0.026) at 90% confidence level. 
If we fix 4> s = 0, the range \${r) S )\ < (0.011 - 0.023) 
is obtained. The parameterizations of F and Fs used 
by the CLEO collaboration 0] yield a comparable limit 
1^(775)1 < 0.013. These results are about one order of 
magnitude more restrictive than the previous best upper 
limit J Sj( 77s) I < 0.19, obtained by the CLEO collabora- 
tion Q. 

Theoretical predictions for ^(j]s) can be given in con- 
text of a MHDM with three or more Higgs doublets [3,[H. 
In such models 77s is given by [I2J 



Vs 



?n T m< 
^11 



X*Z 



(10) 



if numerically small terms proportional to m u are ig- 
nored. Here, M H ± is the mass of the lightest charged 
Higgs boson and the complex constants Z and X de- 
scribe the coupling of the Higgs boson to the r and v T 
and the u and s quarks, respectively (see [l2[ ). The 
limit |3(?7s)| < 0.026 is therefore equivalent to 



\$S(XZ*)\ < 0.15 



M 2 H± 



1 GeV7c 4 



(11) 



In summary, we have searched for CP violation in 
■ K s ^v T decays, analyzing the decay angular 



distributions. No significant CP asymmetry has been 
observed. Upper limits for the CP violation param- 
eter 3(t7s) at 90% confidence level are in the range 
I^K^/s)! < 0.026 or better, depending on the parame- 
terization used to describe the hadronic form factors and 
improve upon previous limits by one order of magnitude. 
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DERIVATION OF CORRECTION FACTORS FOR 
EXPERIMENTAL ASYMMETRIES 

Possible sources of artificial CP violation (CPV) are 
the forward-backward (FB) asymmetry in the e + e~ — > 
t + t~ production and the asymmetry for the detec- 
tion. For both of these experimental asymmetries, cor- 
rection factors have been prepared using data. 



Forward-backward asymmetry 

The FB asymmetry is measured using t + t~ events 
with t~^~ — > i' t tt ± tt + 7t^ on one side and a one prong de- 
cay on the other side. Using the direction of 7r ± 7r + 7r _ 
system as an approximation of the direction of the r lep- 
ton, the measured FB asymmetry is found to be a few 
percent as shown in Fig. Q] and well described by Monte 
Carlo. In order to take this FB asymmetry into account, 
a correction weight factor is prepared as a function of the 
momentum and polar angle of the 7r ± 7r + 7r _ system and 
is applied for each event. The net effect of this correction 
on the measured CP asymmetry is very small [O(10 -4 )]. 
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FIG. 1. Forward-backward asymmetry as a function of the 
polar angle of the three pions, 63^, in the CM frame (z-axis 
along the direction opposite to the positron beam) for t + t~ 
events with t — > u t tt ± ty + tt~ on one side and a one prong 



decay on the other side, 
and Monte Carlo (MC). 



The asymmetry is shown for data 



This detector induced asymmetry is measured using the 
same data sample as for the FB asymmetry studies. The 
charge asymmetry 



A, 



dot l 



^7r : Vir ) 



is measured as a function of the polar angle 9 n and the 
momentum p n of in the laboratory system, where 
denotes the number of particles in each 9^ and p v 
bin. The observed detection asymmetry is up to a 
few percent as shown in Fig. [5] and taken into account by 
applying correction weight factors 

Wdet( fl *>P*) = 2n ± 

for each events as a function of the charge, the mo- 
mentum and the polar angle of the in the signal 



events. The net effect of this correc- 



tion on the measured CP asymmetry is of the order of 
[O(10~ 3 )] as shown in Table I of this publication. 
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FIG. 2. The 7T detection asymmetry as a function of the 
pion momentum ln(p^ b c/l GeV) in the laboratory frame for 
data and Monte Carlo (MC). The asymmetry is shown after 
applying the correction for the forward-backward asymmetry. 
A logarithmic scale is chosen for the momentum in order to 
obtain a more even distribution of events. The current MC 
gives a poor description of this asymmetry because the nuclear 
cross sections for are not well known. 



CP ASYMMETRY IN THE CONTROL SAMPLE 



Asymmetry for the detection 

The different nuclear interaction cross section for tt + 
and 7T~ can lead to a difference in tracking and parti- 
cle identification efficiencies depending on the polar an- 
gle and momentum of the pion in the laboratory frame. 



In order to check for any remaining detector introduced 
bias after the corrections for the FB asymmetry and the 
tt^~ detection asymmetry, the CP asymmetry is measured 
using the data in the control sample. The measured CP 
asymmetry in this sample before and after applying the 
corrections weights as well as the number of the events 
in the control sample are given inTable H 



2 



TABLE I. CP asymmetry A cp measured in the control sam- 
ple for four bins of the hadronic mass W. The 2nd and 3rd 
column show the observed asymmetry with statistical errors 
only, before and after correcting for the FB asymmetry and 
tv^~ detection asymmetry effects. The 5th column shows the 
number of selected events in the control sample. 



W 

(GeV/c 2 ) 



A cp (1(T 3 ) 
Observed Corrected 



Number of events 

m 



0.625-0.890 
0.890-1.110 
1.110-1.420 
1.420-1.775 



-7.3 ±2.0 -1.8 ±2.0 

-4.4 ±0.9 -1.0 ±0.9 

-0.3 ±1.0 1.1 ±1.0 

3.2 ±3.9 3.0 ±3.9 



129993 ± 361 
529559 ± 728 
367181 ± 606 
20274 ± 142 



The resulting C(Q 2 ) can be parameterized as a 7th- 

TABLE II. Parameters describing the parameterization of the 
function C. The coefficients a,i are obtained by fitting C with 
an 7th-order polynomial. The coefficients are obtained from 
the error matrix of the fit. 
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PARAMETRIZATION OF C(Q 2 ) IN EQ. (9) 

The measured CPV asymmetry A cp can be related to 
the CPV parameter 3(r?s) by 



(cos /3 cos V^- — (cos /3 cos + (1) 

COS Pk COS tpk 



— ^'cos&cos^ - 

1 k£T + 



= — I cos Pj cos V*j — cos flj cos "0j 

\jer- jer+ 

—^-T f[[ Q2 ' <Q 2 , cos cosO) 
rii e to d J J Jqi . 

"dr(r-) dT(T+) 



Hi N s 



x cos (3 cos tp 



dio 



div 

C(Q 2 )^^dQ 2 , 



(2) 



with doj = dQ 2 dcos8dcos (3. Here n.; — (n~ + nf)/2, 
where nf denotes the observed number of — > K^-n^Vr 
events in Q 2 bin i (Q 2 e [Q 2 ^, Q%,i\) and N s = J2i n i- 
The coefficients e tot and e(Q 2 , cos/3, cos ip) are the to- 
tal and the three dimensional detector efficiencies and 
r is the total r ± — » J^tt^i- decay width. The func- 
tion C(Q 2 ) contains the model independent terms and 
detector efficiency effects and is obtained after numerical 
integration over cos (3 and cos 9: 



C (Q 2 ) = _I^L s i n 2 



1 
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2\2 



T2m T ~~" ^(4rr) 3 Q 2 
e(Q 2 , cos (3, cos-0) 
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(3) 



etot 



cos (3 cos tp rfcos ^ cos 



For a fixed value of Q 2 , the relation between cos# and 
cos ip can be obtained from Eq. (2) and Eq. (5) of this 
publication. 



order polynomial, 



=,2; 



(4) 



where Q 2 denotes the dimensionless value of Q 2 measured 
in units of GeV 2 /c 4 . The obtained fit result is shown in 
Fig ©. The goodness of the fit is x 2 /ndf = 101.4/92 
where ndf denotes the number of degrees of freedom. 
The uncertainty ac introduced by limited Monte Carlo 
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FIG. 3. Fit of the function C which is used to related the 
measured CPV asymmetry A cp to the CP violation param- 
eter Q(ris)- The function is obtained after numerical inte- 
gration over the cosine of the decay angles (3 and 9 and in- 
cludes detector efficiency effects. The goodness of the fit is 
X 2 /ndf = 101.4/92 where ndf denotes the number of degrees 
of freedom. 



statistics and binning effects in the determination of the 
detector efficiencies can be given as a function of Q 2 



, c (q 2 )=|EEqX 

i=Q j=0 



cov n 2 i 




(5) 



which is obtained from the covariance matrix M?° v of the 

"j. . 

fit. The coefficients ai and ei are given in Table HD 



